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Abstract 

We calculate the one loop amplitudes for the two gluon scattering process in light cone gauge with 
fermions and scalars circulating in the loop. This extends the work of [1,2], in which only the gluon 
circulates the loop. By putting all fields in the adjoint representation with Nf — 2, N a — 6, the scattering 
amplitude of gluon by gluon in the special case of TV = 4 Super Yang-Mills can be obtained. The massive 
fermion and scalar with arbitrary representations are also considered. 
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1 Introduction 



In the previous work of [1, 2], we calculated the one loop gluon scattering amplitude in the light cone world 
sheet setup. The main ingredients of light cone world sheet are explained in [1] and the references therein. 
This setup shows us how to organize the calculation in order to cancel the infamous l/q + (gauge artificial) 
divergences in the light cone gauge propagator. In the cases we have calculated, we were able to show that 
the artificial divergences actually cancel in physical processes, hence any prescription(e.g. principle value, 
Mandelstam-Leibbrandt prescription [6]) should work. For the helicity conserving amplitudes, infrared 
divergences were regulated by cutting off l/q + at q + = 0. Not withstanding this non-covariant regulation, 
we still got a covariant probability by including soft Bremsstrahlung and collinear cmission(absorption). 

This paper extends the work of [1,2]: in addition to the gauge field itself, scalar and fermion fields are put 
into the loop. This paper is organized as follows: in section 2, we review the derivation of light cone gauge 
Feynman rules suitable for massless fields and the light cone world sheet. In section 3, the calculational 
techniques of [2] are briefly reviewed. In section 4, the results are listed without the calculational details, 
since most of calculation is done by computers. When we specialise to the case N g = 1, Nf = 2 (Dirac 
fermions) and N s — 6, we find that the gluon scattering amplitude in TV = 4 SYM takes on a very simple 
form. 

In the last part of the paper(section 5 and 6), we deal with massive matter fields. Some of the tech- 
niques that are designed for massless fields in the first part of the calculation are not useful anymore. So 
we simply follow the textbook procedure, but spinors in light cone are extensively used to organize the 
calculation(section 5). Section 6 lists the results. 



2 The Light Cone Setup 

The light cone gauge Feynman rules are usually obtained by the lagrangian method, namely, we first set 
A- = and integrate out A~ . Then the Feynman rules can be read off from the lagrangian, which is a 
function of A 1 and A 2 . 

Here, a more flexible method is used. The gluon propagator in light cone gauge is: 



fc+ 

k 2 + ie 

Note our metric is diag(l, —1, —1, —1). The numerator can be factored into: 



(1) 



- ^ - ^) = -(4ft + ejtf) - 9 -^k 2 (2) 
Where e AjV are light cone gauge polarization vectors given by 

p _ 1 ^ - ip 2 _. p 1 - ip 2 _ 1 p A _. _p A _ t 

y/2 p° + p A p° + p A V2 P + P + 

They satisfy e • e* = — 1 and k ■ e = 0. These polarisation vectors are defined both on-shell and off-shell. 

The Feynman rules are obtained by dotting the polarisation vectors into the covariant three or four point 
vertices. The tri-gluon vertex Fig.l, for example, becomes 

V ggg = 9f abc (-e^)(-e* 2 „)(-e* 3p )[g^( Pl - p 2 )" + g" p ( P 2 - PzY +9 P «( P 3 - Pif] 

Setting ei, e 2 — e A , e 3 = e v , the above becomes: 



9f abC [(P2 - P3 ) + 4 - (P2 -P3) A + (P3 -Pl) + 4 - (P3 -Pi)"] = 29f abc {Pl t P + + K 2l (4) 
Pi Pi PlP2 
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Figure 1: 



where K?j := (p^Pj — P^pf)- They are related to spinor products according to: 

K i,3 =PiPj\Pi\Pj} =PiPjPiPja 

The spinor notation here is also different from the conventional one [4] . The reader can refer to the appendix 
for an explanation of the spinor notation. 

The gluon propagator Eq.l almost factorises into the product of two polarisation vectors. While the third 
term on the rhs of Eq.2 will make an extra contribution to the four point vertex. For example, consider the 
t-channel diagram Fig. 2, the first two terms of Eq.2 can be associated with the two tri-gluon vertices, the 
third term, which describes the mediation of A~ , gives: 




Hi Pi v, p 2 



Figure 2: 



* * v<nia i * vvp0 5 ttf(pi+Pi) 2 _ iel ■ eXjpj - pf)e* 2 ■ e* 3 (p+ - p+) 

(Pl+Pi) A * {Pi+Pir (Pl+Pj) 

What's happening here is that the explicit factor of k 2 in the third term of Eq.2 cancels the propagator, 
effectively making a four point contact vertex. 




The fcrmion-gluon vertex is the usual igj^t 11 . We can set /i = A or V by dotting — e* v , into /i. 2 

Calling polarisations by V or A is potentially confusing, especially if you arc looking at the diagram up-side-down. 
So, sometimes it is clearer to associate A with 'in' and V with 'out'. 
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Multiplying the spinors to the gamma matrix, we get (assuming now the fermion is left-handed): 



ig(t a ) bc j2ptpi(-eT) [0 pf ] 











Pi 




= -tg(t a ) bc \/2 P +p+{-V2)[ P2 \v](q\pi) 

(6) 



,a Pt 



In Eq.6, in order to avoid defining what is \Jp~^, we choose to associate p\ instead of \jv\Pt to a vertex. 
This won't cause any problem, since fermion loop always closes, so the square root always appears in pairs. 
The fermion propagator is given by: ip^f^ / (p 2 + ie). We can decompose p ■ 7 according to 



>■ a = V2 



p ■ a 



= V2 



-p v p+ 



P v P~ 



= V2p+ 
V2p+ 



-E 







9. 





V2p+ 



V2p+\p](p\ 



V2p- 



-M(v\ (7) 



So the fermion propagator almost factorises too. It is also possible to contract a pair of vertices here: 
the second term of Eq.7 will again cancel a propagator. For example: 



CT, V, Pi 




P, A, p 3 




b, Pi 



c, Pi 



-9\t\ e (n eb fipM\p2\(-e* 3 ■ a) i{P { ^ + + P ^ )2 a (-el ■ a)\ Pl ) {^p\ 

- g \t% e (e\ b ^p~tp~i\p 2 \(V2\ v }{p 3 1) • - ( i V2 • ( t + P * )2 + A v)iv\) ■ (V2\ Pi ]( V \) \pi) 
v (pi+p^r vAPi+pi) 

Pi Pi -+ 



= -2ig\t d ) ce (t a ) eb \ PlP2 + -> -2ig 2 {t d ) ce (t a ) eb ^^ (<S) 
Pi +Pj Pi+ Pi 

The scalar Feynman rules have no suspense in them at all, and can be read off from any field theory 
book. The Feynman rules that pertain to our calculation will be summarised in the appendix. The main 
property of the Feynman rules above is the absence of p~ . 



3 Brief Description of Calculational Procedure (massless case) 

In the above figure 3, the fc's and q are the dual momenta. They are related to the real momenta coming into 
the three legs according to: pi = ki — ko, p 2 = k 2 — ki and p$ = fco — k 2 . The unregulated integrands have a 
symmetry under fcj — > ki+a, which would ensure that each diagram depends only on the real momenta. But 
here, as in [1,2], we use a regulator exp {— ^[(g 1 ) 2 + (q 2 ) 2 }} that breaks this symmetry. Hence, a regulated 
amplitude can depend on the individual dual momenta. This seemingly unwieldy regulator is designed for 
the world sheet purpose, which doesn't concern us here. The result doesn't differ too much when a cut off 
regulator is used. 

The calculation roughly goes as follows, 
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Figure 3: Here we assume only (ko — fo) 2 / 0, and k£ > k± > 



1. Exponentiate all propagators according to 



p 2 + ie 

For the schematic diagram Fig. 3. We have 



dTe iTp2 



~ J W 



4 q i 



) 4 (q-h) 2 (9-fci) 2 (q-k 2 ) 2 
9 dTudTidTtd To{q - ko)2+tTliq - kl)2+lT2{q - k2)2 



(27T) 



2. Integrate out q , leaving a delta function relating g+ to the Feynman parameters (this step requires 
the absence of q~ in all Feynman rules). 

D iT ( 9 -feo) 2 +iT 1 (9-fei) 2 +i7 , 2(g-fe 2 ) 2 



3. Integrate q , q 2 using exp [— <5q 2 ] as a damping factor. 



J dT ° 



Here is the rub: we cannot simply integrate over Tq,T\ and T2, because the prefactor of this diagram will 
have up to second order poles at q + = kf , i = 0, 1, 2 . In order to show the cancellation of gauge artificial 
divergences, we proceed as follows: first eliminate one Feynman parameter in favour of q + : 

k+ < a+ <k+ ■ Ti - ~ fc " +) T - 0- ~ x )b + ~ k o) T _ 1 T T 
fc <q <fc x . J1--7+ p— , J 2- 7+ 7+ . J0-I-J1-J2 

ft l ""2 "'0 

, , . (l-x)(«+ -fc+) ar(g+-fc+) 

fc+< g +<fc+: T 1= ^ f T = ^_ r , IW-Tr-To 

ft l ft 2 ^0 ""2 

Now dT dTidT25(J2Ti — l)<5(X)Tifc i f — <j + ) = dq + dx.J. After integrating out x, we are left with a function 
of <j+ which is defined differently in different regions: k$ < q + < k^ and fc+ < q + < k£ ■ All these can be 
visualised very clearly when we represent a Feynman diagram on the light cone world sheet. The details can 
be found in [1,2]. 

Our observation is that: in each region, all poles cancel. 3 Hence we can perform the final q + integral and 
obtain the results. 



3 in the case helicity conserving amplitude, all poles cancel up to infrared terms, but since infrared divergence is always 
proportional to a tree, they are easy to recognize and deal with. 



4 



4 Results 

Since the Fcynman rules of fermions and scalars are similar to the pure gluon case, the evaluation of scalar 
and fcrmion loop is similar to [2]. In fact, for the helicity conserving amplitudes, the reader needs only work 
out the prefactor of each diagram (as a rational function of q + ) and read off the results from [2] . All of the 
calculation is done by MATLAB program, so we only list the results here. 
We will be following [3] and decompose an n-particle amplitude into: 

M„= Yl Tr(t a H a \..t a ")M(p 1 ,e 1 ;p 2 ,e 2 ;...;p n ,e n ) (9) 

perm' 

where perm' is over non cyclic permutations for complex representation, non cyclic and non reflexive per- 
mutations for real representations. 

4.1 Self-energy Diagrams 



p ki 




Figure 4: A helicity violating self-mass diagram 
The helicity violating self-energy diagram is given by: 

nr = ^^Tr^t 6 ] + fc x + k f] 

n « A = li& TT[tatb] [ fc ° A2 + w + k ^ (10) 

This diagram describes the amplitude of a left handed gluon flip into a right handed gluon. It is only nonzero 
because the regulator we had used doesn't respect Lorentz invariance. This is purely an artifact, and it has 
to be cancelled by a counter term. Also, we observe that ITg A x N s + Up A x Nf + Uq A x N g = in the 
N = 4 SYM case. 

p h 

A, a ^ V,b 
Figure 5: A helicity conserving self-mass diagram 



ttAV 



16tt 2 



a tb~\ 



Tr[t a t 



dxx(l — x)p 2 logx(l — x)p 2 Se r 



-2^TrM 



7 + / dx \x 2 + (1 - xf] p 2 logx(l - x)p 2 5e 1 
o Jo 



dx 



1-x 



1-x 



+ x(l - x) 



p 2 logx(l — x)p 2 8e 1 



(11) 
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In the above list, we have omitted part of the gluon energy terms that is non-covariant (proportional to p + ), 
which was interpreted as a world sheet boundary cosmological constant in [1] . 4 

When we try to fit the self-energy diagrams into the big picture as in Fig. 6, the trace factor becomes 



a d 




b c 



Figure 6: 



(assuming that t a is in a real representation): 

f abe Tr[t e t f ]f fcd = Tr[(-i)[t a ,t b ]{-i)[t c ,t d ]\ = -Tr[t a t b t c t d ] + Tr[t a t b t d t c ] + Tr[t b t a t c t d ] - Tr[t b t a t d t c ] 
-» -2Tr[t a t b t c t d ] (12) 

4.2 Triangle Diagrams 



V, c, p 3 




A, a, pi A, 6, p 2 

Figure 7: Triangle diagram 



The result for the triangle diagram Fig. 7 is: 



pAAV 

1 s 



pAAV 
1 G 



^-Tr{t a t b t c ] 
8ir z 



Pi Pi 



1 1 2r 7 1 l P\Pt 

logp>T---a- - 



6 



P?> 



+2 



2.1 



'.9 



Tr[t a n c ] 



2 , 2 r 7 , 16 1 ptP2 

- 3 lo gPo Se^ + - - a--^ 



3 P 3 +2 



+ infrared terms 



(13) 



where a = 1 if leg 3 is off shell and otherwise, and p is the off shell momentum. The so called infrared 
terms present in the gluon triangle diagram will eventually be combined with infrared terms from other 
diagrams to become proportional to a tree. They are given in [2] and are not repeated here. 
All of the above also contain an anomalous term: 



T ano = 1 {k A + fc A + fc A } 

4 In the light cone setup, certain tadpole diagrams are ill defined, and are discarded from the beginning. This is the reason 
of the non-covariance in the result. 
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T ano = _^9r Tl[t a t b t c ] h k A + fc A + fc A) 
47T O 

r G ° = ^At a t b t C ]\(^ + + k£) 

here the fc's are dual momenta, they arise because 5 is an exponential damping factor of the transverse dual 
momenta. Had we used a cut off regulator ((q — ki) 2 < A for example) instead of 5, the anomalous terms 
would change according to:(k + ki + k 2 ) — > (k — ki) + (k 2 — ki). These polynomial terms must be 
cancelled by counter terms. 

The MHV triangles (with all three legs having the same helicity) give: 



-pAAA 
1 S 



pAAA 
1 F 



-pAAA 

1 a 



87r P1P2P3 

ig Tr\t a t b t°] 

4tt 2 pfPiPt 

^9 rpr.q.b.cl i K 2l) 

s- 2 [ l P tPiPi 



3P 2 



Zpl 



3P 2 



(14) 



4.3 Scattering Amplitude 

Finally we simply list the results for the gluon scattering amplitude with massless fermions and scalars. 
The tree level amplitude for four gluons with the same helicity is zero. At one loop, the amplitude is 

^AAAA = W Tl[t a t b t c td] ^43^32^21^14 

6tt 2 PiP2PiPt st 

_9,Vj 4 ISA ISA rsA ISA 

^aaaa _ ffg Tr[t a t b t c t d ] 43 32 21 14 

37T 2 pi pi pi Pt St 



.4 

-Trr^t c t d l- 

Pi Pi Pi Pi st 

The tree level amplitude with only one unlike helicity is zero too. At one loop level, the amplitude is 



j n 1 If A ISA tsA tsA 

^AAAA = 19 Tl[t a t b t c td] ^43^32^21^14 ^ 

3n vtvtvtvlst 



l9 \TT[t a t b t c t d }(s + t)- K ™ p Z p * 



487T 2 L JV 

A AA AV _ l 9 T r.a.b.c f dy , ,N 13^2^4 

" ^ ' * ] (S + ^ ^3^21*14 

^ AAV - ^Trr^]( S + ^C P lvV < 16 > 

Z47T /Y 43 l\ 32 l\ 21 _f\ 14 



The helicity conserving amplitude is non-zero at tree level. They are given by [4]: 

_o J^A4 + + 

MAW _ - 2 fabe feed zrL 12 Ps Pa 

-^tree L 9 J J ry/\ is A is A ts A , 



K^K^K^ptpi 
aAVAV - 2 fabe feed - 2K uPipi 

5/ ' ^^^i^Vi^ { } 

The factor f abe f ecd can be converted to -l/G(G)Tr , i 6 ] [i c , -» -2/G(G)Tr [t a t4 c i d ] . 
At one loop level, the amputated Green's function is (with infrared terms omitted) 



87T 2 H#&# 3 VMipM 



^ + ilog^ 



ix + I 

6 3 



7 



A AAVV 



A AAVV 
A G 



An 2 



Tr[t a t b t c t d ] 



19 2 1 X 7* 



X 



^ Tr[i a i b £ c t d 1 < ^^12^3^4 

8tT 2 J \if^JfAjfAp+ ft + 



(log : 



£ + 7T 2 ) - ill og^i+ 7 4 
t ' 3 9 



(18) 



4 AVAV 



4 AVAV 



A AVAV 



^g 4 Trff a f b i c i d l I ^-^"i3 4 P2"P4 

87T 2 I /T A /T A /T A ?T A n + n + 

on lJ\ i3 l\ 32 J\ 21 J\ 14 p 1 p 3 



sH 2 



+ 



s(2t 2 - 5st - s 2 ) 



-ig_ 
An 2 



6(s + t) 3 



Tr [t a t h t c t d ] 



log : 



^ l0g,5e7s+ 2(7T^ 



2(s + t) 4 
ts 1 
18 



(lQg^+TT 2 ) 



t 



6 3 



+ 



s(5t 2 + si + 2s 2 ) 
3(s + t) 3 



log - — - log(5e 7 s 4 
6 t 3 6 



2(s + t)* (lQS t 



ts 



(s + i) 2 



19 
~9~ 



1} 



■ 4 

8tt 



-2K^+ 



13P2P4 



+ 



s (Ut 2 + 19st+ lis 2 ) , s 11 
: loe 



(s2 + St + t2) -(log 2 £ + . 2 ) 



(s + i) 



■ log <Je~ 



(t + s)4 

ts 



73 



1 



(19) 



The symbol X above is the relevant four point vertex: — 2(p J [p 3 + P2pt)/[(pt + Pt)(?2 + P3)] or 

2(ptP3 +PiPi)[(Pi +Pt){pt +Pa)\ + 2 (PiP2 +P3P4 )/[(Pi" +Pa ){P2 +Pt)\- Thc four P oint vcrtex is not 
a valid counter term, because it is not a polynomial of momenta. But we can add a term proportional to p 2 
to the self-energy term Eq.ll, this only changes the field strength renormalisation by a constant, hence is 
perfectly allowed. With this term, the coefficient of s and t channel exchange diagram is shifted. So if we 
pick the numerical factor in front of p 2 to be — |, — g, — i for scalar, fcrmion and gluon respectively, then 
they will match the coefficient of the lone four point vertex, completing it to a full tree. This brings about 
a change in the numerical factor:-^ — ► — yg, ^ — > ip — > ^ [2]. In fact, thc numerical factor is quite 
unimportant, as it can be altered by a redefinition of coupling constant. 

It is tempting to go to TV = 4 SYM by putting N s 
is amusing(up to infrared terms): 



6, Nf = 2(Dirac fermions) and N g = 1. The result 



/l A AAA 
A SYM 

A AAVV 
A SYM 



-1 AVAV 
A SYM 



/l A AAV _ n 



SYM 
4 



Tr[t a t b t c t d ^ 



8tt 2 
• 4 

^Tv[t a t b t c t d \ 
8n 2 K&K&K&Khp+Pt 



-2K^p+p+ 



-(l0g 2 £+7T 2 ) 

-(log 2 ^ 2 ) 



(20) 



Quite remarkably, in this case, all counter terms cancel. 



4.4 A Word on Infrared Terms 



In the above list of results, we have omitted the infrared sensitive terms. First, as we can see that the 
self-energy diagram always contains a term log(5e 7 p 2 x(l — x). This gives a multi-particle branch cut on the 
positive real axis, and stops us from doing wave function renormalisation. This can be cured by summing 
over collinear emissions or absorptions. The analysis of [2] showed that doing so is equivalent to replacing 
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log<5e 7 p 2 x(l — x) with log(5e 7 A 2 x(l — x), A being the jet resolution. For fermion loop or scalar loop, this 
substitution alone is enough to regulate the infrared divergence (the triangle or box diagrams involving 
fermions or scalars are devoid of further infrared divergences). While for gluons, we refer to [2] for a 
complete treatment of infrared terms. Here we only list the infrared terms from triangle and box diagrams 
from the gluon loop. The k below is dual momentum, with p^ = k^ — k$, p\ = k^ — k^ , p 3 = k 3 — k^ , 



pi 



fcg" — fcj" and fcg" < k 3 < k^ < k^ ■ Pi P2 are incoming legs while pz P4 are outgoing. 



k+ <q+ < k+ 



k^ <q + < k^ 



kf <q + < k^ 



+ 
+ 
+ 



+ 
+ 



+ 



+ 
+ 
+ 



^L-Tr[t a t b t c t d ] x 

871" 



Og 
Og 

og 
og 

og 
og 
og 
og 



(fc+ - q + )(-k+ + q+) 3 6e> 

(k+-k+r 

(k+ - q+)(q+ - fc+)sfe 7 

(ki — k^){k2 — k$) 
(k+ - q+)(-k+ + q+) 3 6e* 

(-k+ + q+) 2 t8e> 



1 



r - k 3 

2 



-kf 



q+ - k+ 



(k l k )(k 3 k ) 

(k+ - q+)(-k+ + q+)sSeJ 

(k+ k+) 2 
(k+ - q+)(q+ - kt)t8e> 

{-kt + k+f 

(k+ - q+)(q+ - k+)s5e"< 

(—k^ + k^)(kQ — k^) 
(g + - k+)(k+ - q+)aSe-y 

(^2 ~ ^o)(^t — ^3~) 
(-k+ + q+)(k+-q+)t5ei 
(k 3 — fc^)(— fc^ + fcj") 



q+ - fc+ 



k+ + q 



k,2 



q+-k+ 



k 3 



-){q+ - k+)t5e~< 



k+ 



(-k+ + k+)(-k+ + k+) 

(k+ -q+)(-k+ + g+)sfe 7 

(k^ — k^) 2 
(q+ - k+){-k+ + q+)s5e~t 

(k+ -k+)(k+ -k+) [ q + - /,; 

(q+ - k+)(k+ - q+)s5e~< 
(^2 ~ ^o)(^t — ^3~) 

(-fc+ +g+)HSe^ 
{k,2 — k^)(k2 — k 3 ) 



k 3 



-*+-k+ 



k$ 



The above is multiplied by the corresponding tree amplitude (either — 2Ki 2 4 piP4 / {K^K^K^K^p^p^) or 

-20+P4 + /(^4 A 3^ A 2^2i^i>|P^) A A VV and A V AV). 

The infrared terms above will be combined with soft Bremsstrahlung and collinear emission (absorptions) 
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along with the self-energy insertions on the external legs to give a finite result. Note also that the log(5e 7 
part of the above infrared terms actually cancel the divergent x integral in the log Se 1 part of the the gluon 
self-energy diagrams on the external legs. Thus the coefficients of logfc 7 become 1/6, —2/3 and — 11/3 (times 
tree) for scalar fermion and scalar respectively. Here we quote the scattering probability with general massless 
gauge and matter fields, with initial and final states treated as in [2]. 



\A A 



1 + 



g 2 c(G) 

8tt 2 



„, 2 A 2 „, o A 2 2tt 2 r 67 5 Ar 26 Ar . 

2l0g v- 21 ° g V + — + It-18^"T^ 



, (5e 7 A 4 r 11 1 T 4 Ar . , 2 s' 

+ log T [ -y + 6^ + 3^ ]+los iti 



l AVAVI2 



1 | a 2 c(G) 



log 



8tt 2 



5e~*A 4 r 11 1 



„, 2 A 2 „, 2 A 2 2tt 2 r 67 5 Ar 26 Ar . 

2l0g v- 21 ° g V + — + It-I8^-t^ 



1 



, - [ -y + 6 Ars + 3^ ] + 2 (s+ t )4 - \t 



log 2 — [s 2 i 2 iV s + 2st(s 2 + t 2 )N f + 2(a 2 + si + i 2 ) 2 ] 



1 



+ — -3 log r A-s(2t 2 - 5st - s 2 )N s + 4s(5t 2 + st + 2a 2 )N f - 2s{Ut 2 + 19 st + lis 2 )} 

6(s + t) \t\ 

st 



2(s + t) 2 
where A is the jet resolution 



N s + 4N f - 2} 



(21) 



5 Massive Matter Fields 

The main difficulty for massive field calculation is the Feynman parameter integrals. Indeed, we can only 
reduce all Feynman parameter integrals into a set of three definitive integrals. They are: 

1(a) 



So dX sx(l - 

f 1 1 
J(s) := / dx— log 

Jo x 



-x)+M 
sx(l -x) + M 
M 



K(s t) ■ I" dx = log isxjl - x) + M)(tx{l - x) + M) 

[ , ) - J aX stx(l-x)+M(s + t) ° g M 2 {ZZ) 

Where M is in fact -to 2 + it. As M — » 0: 

1(a) ~ 21og^ 

J(a) ~ ^log 2 ^r 
w 2 M 

K(s,t) ~ -. 2 + 21og^log^ 
K(s,t)-2J(s)-2J(t) ~ -7r 2 -log 2 ^ (23) 
For example, a simple integral: 



I 



5(Ti +T + Tz + T 2 - VdT^TodTadT: 



TiT (l - Ti - Tp) 

Iq „ v _ , _ , _ , _ a _ 2 (TlT3t + ToT2S + M) 2 

M M M 1 

can be reduced to: — ^ r J(s) tt, r J(t) — - — 7 -K(s,t) + - — (24) 

s 2 (s + t) w t 2 (s + t) w 2st(s + t) v ' ; 2st K ' 
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The Fcynman rules for scalar fields remain good, but the decomposition Eq.7 does not, since it introduces 
l/<7 + factors into the Feynman rules. This will complicate the already complicated Feynman parameter 
integral. Also q + and q~ will be treated the same in contrast to the massless case, where q~ is integrated 
out and q + is given by ^2 Tik^~ / ^2 ^i- 

In order to organize the gamma matrix algebra in the fermion part of the calculation, we make use of 
the factorisability of gluon polarisation vectors Eq.38 to reduce products of gamma matrices to products 
of Kjfj — p^Pj — Pj~Pi, which had been proved to be quite handy. For example, if we are to calculate the 
diagram Fig. 4. We would write down: 



(tgYTT(t a t b )TT ([(q - fc ) • 7 + m]^[{q - fci) • 7 + m]>) (-#)(- e^) 



(q — fco) 2 — m 2 (q — fci) 2 



(25) 



The numerator can be written: 
Tr 



m (q — fc ) ■ a 

(q — fco) ' v m 

m (q — ki) ■ a 

{q — k\)-a m 



-V2\v)[k! -k \ 

-y/2\ko- 

-x/2|fc — fci>[r?| 

-v^N<fci-*o| o 



= [v\Q - ko\v)[ki - k \q - fei| A; - fci) + (fco - h\q - fc |fci - fc ] (r)\q - fci|?j] 
+ (fc - h\m\r]} [fci - k Q \m\rj\ + [r?|m|fci - k ](r]\m\k - fci) 



(26) 



For the notation of spinor, see the appendix. The standard procedure of momentum integration tell us to 
shift q — > q + xk + (1 - x)k\, and to replace q^q" by q 2 g^ /4, etc. 5 : 



[r)\q - fco|»7>[fci - k \q - fa\k - fci) 

= [?j|(l - a;) (fci - fco)|»?)[fci - fcoKfco - fci)|fco - fci) + [»?|^ a |j?>[fci - fcol^Jfeo - fci) 
= (1 - x)V2{h fc ) + ^{k>-kl)+ + Mh k ° ]{k ° ~ kllv) T 

2 



Ab 



= —x(l - a;) (fco - fci) 2 
here we have used an identity: 

VaaPhbV^ = 2e db e *b, CT M (7^, .9^ = 2e £ 

For a more complicated example, a string of spinor products becomes, 

[PiW " k \p 4 )[r]\q - k 3 \p 3 ) 
= [ Pl \q + Ti(fci - fc ) + T 2 (k 2 - fc ) + T 3 (fc 3 - fco)|p 4 ) fr|<Z + ^o(fco - fc 3 ) 



(27) 



(28) 



Ti(fci-fc 3 )+T 2 (fc 2 -fc 3 )b3) 



tPiN<P3|P4>2- + bi|72P2|p4)[?7|7bp4 - 7\p 2 b3) 



1 



V2T 2 

P1P2P4 



-i\ 12 i\ 42 ^ i A. 



P3 + 



34 



+ — ^32 J 

P3 



Basically, all of the spinor products can be reduced to one of the K^s, and the reader can find in the 
appendix some practical details as to how to organize the products of K^s. After the momentum integral 
is done, we can perform the Feynman parameter integrals using Eq.22. 

5 We have omitted the complication due to the <5 regulator. For one thing, the momentum integration is no longer homoge- 
neous, hence the replacement of q^q v by q 2 g^ /4 is problematic. Some extra work is required to sort this out, which we will 
not bore the reader with. 
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6 Results 



6.1 Self-energy Diagrams 

For Fig. 4, the results are: 



2 1 x 7 2 . 1 (P - 4 ™ ) j, 2n 5 2 4 2 

p"log5eW + — — I(jp z ) - —p z + -m z 

2] A7 2 l (p 2 -4m 2 )(p 2 + 2m 2 ) „ 2 , , 13 2 , 8 2 
p logde'm — — l{p) + —p + —m 



3 p 2 9 r 3 

The gluon mass counter term in the expressions above has been removed already. 



(29) 



6.2 Triangle Diagrams 

The results for triangle diagram Fig. 7 are: 



pAAV 

1 s 



-pAAV 
1 F 



07T Z Pi P2 



ilogm 2 ^ + (p 2 -4m 2 )( ry 



(p 2 — 4m 2 ) P1P2 to2 



P3 +2 



K(p2) 



pfpt m 2 . 1 4m 2 Iptpt ,„ 24m 2 , 

— J Po - s + Tj- + —5-) 

P3 9 3p 2 6 p+ 2 p 2 



-T-I^ a t h t c \ x —"21 



£ logm 2 ^ + (p 2 - 4m 2 )(- (P ° + 2 4 m) + 2«^#^)/(p^) 



3P* 



pr 



16 8m 2 lPi"P2" n 24m 2 , 



,+2 „2 



(30) 



where a = 1 if leg 3 is off shell and otherwise, and p a is the off shell momentum. Their anomalous terms 
are identical with the massless result, which agrees with the fact that anomalous terms are UV effects. 
The result for the MHV triangle is: 



-pAAA 

1 s 



pAAA 
J- P 



\+»+hc] ^2T) 3 [ 4m 2 (p 2 - 4m 2 ) J( 2 , , 4m 2 J( ^ , 2(p 2 D + 24m 2 ) 
8tt 2 [ 

4tt 2 



P1P2P3 



pl 



-Hpi) + —J{p z o) + 

Po 



Pi P2 P3 

6.3 Scattering Amplitudes 



8m 2 (pg-4m 2 ) 2 ^ , 8m 2 2 ^ , 4(pg + 24m 2 ) 
~ % KPo) + ^JiPo) + , pt 



(31) 



^AAAA = l 9 T ^ t a t b t c t d^ 
7T 

A^ A ^^Tr[t a t b t c t d ] 

7T 



K A K A K A K A 
^ 43 ^ 32 lx 2\ sx 14 


"1 


4 

m 


pi pi pi pi st 


6 " 


" ~~st 


K A K A K A K A 

- tY 43 /Y 32- tl 21- tl 14 


"1 


2m 


pi pi pi pi st 


3 " 


st 



-K{s,t) 



(32) 



^AAAV 



43- ,1 32 J1 21 

(2i + s)tm 2 4m 2 



2(s + i)s 
(2s + t)t 2 m 2 



)J(s) _(2 S + ^m 2 (1 



2(s + 



T/ , (2t + s)s 2 m 2 . , 



stm 



2(s + t) 2 



4m £ 
(1- 



)/(*) 

2 m 2 (s + t) 
si 



MM) 



12 



2(2s 2 - si + 2t 2 )m 2 + (s + t) 



st 



6 



4 AAAV 



-*9 
4tt 2 



Tr[i a ^t c t d ; 



K y K v K A 



(2t + s)tm 2 _ 4m 2 _ (2s + t)sm 2 _ 4m 2 



+ 



+ 



(s + i)s 
2(2s + t)t 2 m 2 

(s + t) 2 s 
A{2s 2 - st + 2t 2 )w? 



J(s) 



s ' - ' (s + t)t 
2{2t + s)s 2 m 2 



+ 



(s + t) 2 t 

( s + ty 



-j(t) 



stm 

Js~Tt) 



t 



r(l-2 



m 2 (s + 1) 



st 



)K(s,t) 



(33) 



In contrast to the massless case, the external leg factors are included in the results below for the helicity 
conserving amplitude. 



A AAVV 



A AAVV 



Stt 2 K&K&K&Kfo+Pt 



m 2 2m 
^~2s~ ~ ~it 



M3 J 

2m 4 



st 



4m 4 1 . T , . 



-^K(s,t) 



2m 2 



4m z 
~3T 



1 

18 



-\ log<5e 7 ™ 2 - \ 
6 3 



4tt 2 



Ti[t a t b t c t d ' 



'12P3P4 



/ m 2 2m A 



4m 



st 



+ 



8m 4 
Itt 2 " 



+ 



. , . 4m 2 8m 2 19 



2 2 

+ - log(5e 7 m 2 

3 6 3 



3 3 



(34) 



^AVAV 



A AVAV 



^ 4 Tr[t a t b t c t d 1 ^^-ijpjpl 

Stt 2 K&K&K&Khp+Pt 

tm 2 (l7st + s 2 + 4t 2 ) 2tm 4 (2i + 5s) t(5si - 2s 2 + t 2 ) 



6(s + i) 3 s 



3(s + t) 2 S 2 



12(s + t) 3 



sm 2 (17si + t 2 +4s 2 ) 2sm 4 (5i + 2s) s(5st - 2t 2 + s 2 ) 

~H a I,. , I\3l 1 oC„ , I\2l2 1 .\3 )■'(*) 



6(s + t) 3 t 



+ ( 
+ 



2stm 2 

JsTtf 

2m 2 (2t 2 - 



s 2 t 2 



(s + *) 4 
2s 2 + 3si) 



3(s + i) 2 i 2 
)(J(*) + J(t)) + (- 



12(s + t) 3 



2stm 2 



+ 



m 



+ 



sH 2 



(s + t) 3 (s + t) 2 2(s + t) 



:)K(s,t) 



3(s + t)st 

■— log<5e 7 m 2 — - 
6 6 3 



+ 



1 



- -X + 
6 3 



s 2 + t 2 + list 
18(s + t) 2 
1 



4tt 2 



Tr[t a t b t c t d ] 



-2K^p+p+ 



tm 2 {bs 2 + 2t 2 -5st) Atm 4 (2t + 5s) 



3(s + i) 3 s 



3(s + i) 2 s 2 



t(5s 2 + 2t 2 + st) 
6(s + tf 



)I(s) 



tm 2 (5t 2 + 2s 2 - 5st) 4sm 4 (5t + 2s) s(5t 2 + 2s 2 + st) 
+ V o7TT~A3+ ' o7TT~iA272 c? „ 1 773 



3(s + t) 3 i 



3(s + t) 2 t 2 



6(s + t) 3 
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Astm 2 st(s 2 + t 2 ) 



){J{s) + J{t)) + { 



m 2 (s - if 



+ 



2m 4 st(s 2 + 1 2 ) 

(s + 1) 2 ~ 2(3 + ty 



)K(s,t) 



+ 4m 2 (2t 2 + 2.s 2 + 3si) 19s 2 + 19< 2 + 47st 



3(s + t)st 
+ ^log5e 7 m 2 - | 



9(s + t) 2 



1-1 



(35; 



The last numerical factors —1/3 and —2/3 inside each square bracket is due to the amputation of external 
legs: [lim p2 ^ n(p 2 )/p 2 ] 1 /2. 



7 Concluding Remarks 

In summary, we have calculated the next to leading order gluon scattering amplitudes with fermions and 
scalars in the loop, as an extension of [1,2]. They are all part of the project of showing that summing planar 
diagrams on the light cone world sheet does give correct results (at least at one loop order). The light cone 
world sheet tells us to organize diagrams such that an explicit cancellation of gauge artificial divergence can 
be seen. There is no gauge artificial divergences in the fermion and scalar part of the calculation, and the 
method developed in [2] works out as expected. And the same counter terms required in [2] appear here too. 

In the special case of Af = 4 SYM, all of these counter terms cancel among different sectors. This is a 
mere additional bonus from the field theory point of view, but quite crucial for the project proposed in [1, 2]. 
More specifically, in [2] , the dimension of the target space was increased to accommodate extra world sheet 
fields (spurions), these fields couple to gluons in such a way that will produce the counter terms needed (the 
four point contact vertex and a 2/3 in Eq.19). In order to formulate N = 4 SYM on the world sheet we 
need to increase the space time dimension to ten, so as to have six more world sheet fields in addition to the 
original q*(cr, r), i = 1, 2 [5]. These eight fields correspond to the eight bosonic degrees of freedom in TV = 4 
SYM. The AdS/CFT duality says this theory is dual to the string theory in 10 dimensions, meaning we no 
longer have room for extra dimensions to produce counter terms. 

The massive matter contribution to the gluon scattering was also calculated for completeness. An expo- 
nential damping factor was used instead of the dimension reduction as UV regulator, so the box integral is 
slightly different from [7, 8], so no machinery of hypergeometric functions need be used. 
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A Notation 



<j» :=(!,<?) a" :=(!,-<?) 



Pa = e abP b p a = Pb e ba 

\p\~Pa \ P )-=p a [ P \-=p h (p\-=Pa 



?7" = T}" 





-1 



(36) 
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We can define the light cone version spinor as: 



P 



_ P A 




1 


1 


Pa = 


P V 




p+ 



(37) 



The spinors satisfy the Dirac equation if p is light like. Note that they don't have the correct normalisation, 
namely p ■ a aa ^ p a p a , but they have the merit that p a = (—p) a - The polarisation vectors of gluon can be 
written as: 



V2 



\p](v\ 



-V2\p](r)\ 



[v \p] 

= PjPi (Pj \Pi ) , K ij= PiPj \Pi \Pj ] 



V2 



v](p\ 
\p\v] 



-V2\ V }(j>\ 



(38) 



The Kij's satisfy: 



E 



r> + K* 1 
Pi ^jk 



tin 



E 



Pi 



(39) 



j run from 1 to 4, but only two of the six -Ky's are independent, 

4 , 



's with total helicity 4 can be reduced to either (K^ 3 ) 4 or (K^ 3 ) 3 K 32 



In the current case we are dealing with 
say if 43 and K32 . And any product of 
Product of helicity 2 can be reduced to (K^ 3 ) 2 and K^ 3 K 32 . Product of helicity can be reduced to 1 and 
K^ 3 K 32 . The reduction is in general a formidable task for human, but quite a piece of cake for computers, 
as all our calculations are done with computers. 



B Feynman Rules 




The fermion gluon vertex Fig. 8 is given by: 

= -2ig(t a ) bc -^K£ uq V r A = -2ig(t a ) bc ±K^ q 

V? = -2ig(t a ) bc ^K^ q K v = -2i fl (r) fec ^_if p v i>? (40) 

qt- q + Pi 

where the subscript /, r denotes whether the fermion is left or right handed. 
The gluon scalar vertex Fig. 9 is given by: 

= -2ig{t a ) bc ^-^K^2 (41) 
P b +pJ 
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fj,, a, q 




/ \ 
Figure 9: gluon-scalar-scalar 3 point vertex 

a, pi d, p 3 




' Pi c > P2 6, Pl c, p 2 



Figure 10: Two diagrams contribute to the fermion-gluon 4 point vertex. Note these are not exchange 
diagrams but rather four point vertices gotten through cancelling a 1/p 2 factor in a fermion propagator or 
a gluon propagator, see Eq.7 



Here we have used real scalar fields and hence it transforms in a real representation. 
The gluon scalar vertex Fig. 10 is given by: 

{pi - pirn 

n =-Mg (t ) ce (t ) eb — j— — — + zg / [t )cb —— j— + , 2 

Pl + Pi \Pz + Pi I 

t/AV _ r* 2 f dae(j.e\ (pi -pim 

= 29f {t )cb (pi + pi? + 

(pi -pi)pi 
Vl =2gf {t u (Pi +P ir 

where the first super index on V refer to the polarisation of ps, the second p4- 
The gluon scalar four point vertex Fig. 11 is given by: 

V ggss = -z.9 2 [(n ce (t%d + (t b ) C e(t a U + gy abe (nJP^-4^i^ 

(Pl + P2 ) 

For tri-gluon vertex Fig.l 

V ggg = -9f abc [t*i ■ e* 2 (pi - p 2 ) • e* + e* • e*(p 2 - p 3 ) ■ e* + • £l (p 3 - Pi) ■ e£] 
Setting /x, f = A, p = V, the above becomes: 

.9.r bC [(p 2 -p 3 )+4 - -P3) A + (P3 "Pl) + 4 - (P3 -Pl) A l 



A 



pi>2 



1G 



A, b, p4 V, a, p 3 




c, Pi d, p 2 



Figure 11: Scalar-gluon 4 point vertex. Again, the second diagram is not an exchange diagram but rather 
a four point vertex obtained through shrinking a gluon propagator 



P, C, P3 




H, a, pi v, b, p 2 

Figure 12: Tri-gluon vertex 




(j,, a, pi v, b, p 2 M, a , Pi v , b , P2 

Figure 13: Gluon 4 point vertex 



The gluon four point vertex receives contribution from two sources: the left diagram in Fig. 13 is simply 
the covariant four point vertex: 

Vi = *s 2 r b 7 ec W p - 9 w 9 va \ + w 2 f dae r bc [g^g pa - g w g va \ + *g 2 f cae f ebd [g^ g pa - g p °g vp \ (45) 

The second is obtained by shrinking a propagator. 

y 2 = g 2 fd ae f e bc[el , £ , ^ _ ^ *f V + frl ±J*L [ e « . e g - p,) fj] 



■2 fdae f ebc (el \ £ l)( e 2 ' 4) (P4 ~ Pi ) + ~ P 3 ) + 

There are two cases in which the gluon four point vertex is nonzero: 



t 7 _ Q ■ „2 /dae /efec Pl~P3 j" Pi ■ „2 face /ebd P^P^ + gjjj" 

w ~ 5/ 7 (Pt+PiM+Pi) + 19 f f (pi+PiM+Pi) 



17 



,2 rabe feed P2 P3 + Pi Pa 



{pi + pi)(pt + pi) 




(46) 



When using these vertices, we need to watch the indices of structure constants closely: not all terms are 
going to make contributions to Ti[t a t b t c t d ]. 
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